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SUMMARY 


To account for the turbulent teirperature and species-concentration fluctuations , 
a model is presented on the effects of chemical reaction rates in computer analyses 
of turbulent reacting flows. The analysis results in two parameters, called the 
temperature and species amplification ratios, which multiply the terms in the equa- 
tions to describe the rates of production or depletion of the species (i.e. , the 
reaction-rate equations). The amplification ratios are dependent upon the mean val- 
ues and the intensity of the turbulent fluctuations of the temperature and species 
concentrations. Graphs are presented of the amplification ratios which will permit 
use of the model in existing turbulent flow-field analyses. 

The model was used in a two-dimensional parabolic-flow computer code to com- 
pare the results obtained with ejq^erimental data available for an axisymmetric, 
supersonic-combustion, hydrogen jet burning in air. The theoretical results show 
some improvement in agreement with the experiment in the far field (i.e., greater 
than 20 jet diameters from the injection station) in the region near the jet center- 
line. Recommendations for further improvement of the model include additional 
assessment of several of the parameters which arise in the formulation of the model 
by additional comparisons with experimental data. 


INTRODUCTION 

The highly nonlinear coupling of turbulence and chemical reaction increases the 
difficulty involved in the numerical computation of turbulent reacting flows. This 
difficulty occurs when global models of the finite-rate reaction are employed, as 
well as with models incorporating detailed chemical kinetics (ref. 1). The necessity 
for dealing with detailed kinetics when modeling turbulent diffusion flames is made 
clear in reference 2. The problem dealt with herein is the modeling of turbulent 
reaction rates in the numerical computation of combustion flow fields which employ 
detailed chemical kinetic mechanisms. Most such analyses (e.g., ref. 2) employ the 
mean temperature to calculate the Arrhenius reaction-rate coefficient; the coeffi- 
cient is then multiplied by mean values of the species concentrations to produce a 
pseudomean reaction rate. In the formulation developed herein, a turbulent mean 
reaction rate is proposed which involves terms which are the product of a mean rate 
coefficient and the mean of the product of the species fluctuations. Quantities 
termed "amplification ratios" (presented as a series of graphs), which are dependent 
upon the mean temperature, the mean species concentrations, and the intensity of the 
turbulent fluctuations of these quantities, multiply each pseudo reaction-rate term 
and effectively "correct" these terms to their turbulent mean values. The proposed 
model is applicable to global and quasi-global mechanisms (ref. 3), although the 
principal consideration here is the development of a computationally viable approach 
applicable to large chemical mechanisms. 

Among the approaches previously employed is that used in reference 2 in which 
the unmixedness model suggested in reference 4 is applied to terms in the reaction- 
rate ejqpressions which contain either the fuel or the oxidizer concentration. An 
alternative approach (refs. 5 to 7) solves transport equations for joint probability 
density functions (pdf) of the species using Monte Carlo solution methods. Other 
approaches which have been examined include the expansion of the exponential term in 



the Arrhenius ej^ression for the reaction-rate constant (ref. 8), and the asymptotic 
es^ansion approach in reference 9. Heuristic models (refs. 10 and 11) are designed 
to circumvent the requirement for detailed chemical kinetics entirely. 


The model proposed herein is based upon the assumed pdf approach employed in a 
large number of prior studies (e.g. , refs. 1 2 to 1 6) . The particular approach sug- 
gested is described in the next section, where probability density functions for the 
temperature are examined. Also proposed in the next section is an assumed pdf model 
to deal with the turbulent fluctuations in the species concentrations. Results are 
presented using the proposed model to calculate a supersonic, axisymmetric flow field 
for which detailed experimental data were available. 

The contributions of the following colleagues to the work reported herein are 
acknowledged: Paul J. Antaki, Eileen Foy, Deborah Goldstein, Gregory Kassar, and 

Gregory Ronan at The Cooper Union. 


SYMBOLS 


pre exponential constant in rate coefficient equation 
mass fraction of species j ( j = A^ , A 2 , . . . ) 
temperature exponent in rate coefficient equation 


d 

E 


K^,K2,K3,K4 

^b 

kf 

k(x) 

P 

pdf 


Pt 

p(x) 

q 

R 


molar concentration of species j 
diameter, m 

activation energy, kJ-mole“^ 
constants 

backward reaction-rate coefficient 
forward reaction-rate coefficient 

reaction-rate coefficient of fluctuating variable x 

static pressure. Pa 

probability density function 

pitot pressure. Pa 

pdf of fluctuating variable x 

a priori probability in joint most-likely pdf 

gas constant 


2 



r 

T 

t 

w. 

J 

w. 

3 

• 

w. 

3 

X 

Z(x) 

X 

P* 


radius measured from jet centerline, m 
dimensionless molar concentration of species j 
temperature, K 

ratio of activation energy to gas constant, K 
dimensionless teirperature 

instantaneous net rate of reaction of species j, moles-(m) 

— 3 —1 

time-averaged reaction rate of species j, moles-(m) s 

— 3 — 1 

pseudomean reaction rate of species j, moles-(m) s 

distance along flow direction, m 

species amplification ratio 

temperature amplification ratio 

amplification ratio of fluctuating variable x 

coefficient in most-likely pdf 

correlation coefficient 


Subscripts: 


max 


maximum 


min minimum 

n nozzle 

ref reference 

Chemical notation: 




H20 


N- 


hydrogen 

water 

nitrogen 

oxygen 


A prime with a symbol denotes a fluctuating component, 



DESCRIPTION OF REACTION-RATE MODEL 
Assumed Probability-Density-Function Approach 
For the general chemical kinetic mechanism 

(i) A^ + A 2 ^ A^ t A^ 

(ii) A^ + A^ 1: Ag + Ay 


the instantaneous reaction rate of species A^ is 


-JCItv Cjv + k, iCj, Cj. - k^ ..Cj. C, + ku 44C, c, + 

A I/i A.^ A2 bf X A^ A^ f/Xi A-j A^ b/XX Ag Ay 


( 1 ) 


where Cj is the molar concentration of species j# 

In a turbulent flow, the turbulent mean value of w is 


'^A " J J 

^ —00 —00 


/ w p /k . » . 


f • • • / c f c , * * *1 


dT dc 


where •••) is the joint pdf of the indicated variables* Such a formulation 

cannot be modeled at present and is, in any event, mathematically intractable. 

Hence, it is proposed to treat each term in equation (1) separately and then sum, as 
f ollows: 


w, = -w, + w, + . . . 

1 1,1 1,2 


where w = -k c c ,w = -k • ^ / etc., and where, for example, 

1,1 *1 *2 1,2 *3 *4 


w_ 


= / / / w 


1 ^ <-00 -,00 —00 




dT dc dc 

""2 




dT dc, dc, 
""i ""2 


( 2 ) 


If it is further assumed that the turbulent temperature fluctuations are statisti- 
cally independent of the fluctuations of the species concentrations (they are, of 
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course, mutually dependent through the transport equations describing the turbulent 
flow field), then equation (2) may be written 



[/” ^(T) p(T) dTirj” /” < 

00 ' J |— 00 —GO “J 





(3) 


where p(T) is the pdf of the temperature and ^ joint pdf of the 

species concentrations. Although this assumption of statistical independence cannot 
be justified a priori on physical grounds, it is demonstrated later in this section 
that the results produced by a fully dependent statistical model do not differ 
greatly from the model proposed here. On the other hand, the proposed model is, from 
a computational standpoint, substantially simpler to incorporate into existing com- 
puter programs. 

The integrals in equation (3) may be bounded between zero and unity by the 
introduction of the following parameters: The dimensionless temperature 


T - T . 

min 


T - T . 
max min 


(4) 


and the dimensionless molar concentration of species j 


^j ” ^j,min 

c — c 

j,max j,min 


(5) 


where the quantities T^. , T c. and c. represent appropriate mini- 

^ min' max 3 , min j,max ^ 

mum and maximum values of the variables, specified a priori. For example, 
could be the initial reactant temperature and the adiabatic flame tempera- 

ture. Also C;j^niin zero and ^xaax '^ould be the initial value for a reactant 
species, the final equilibrium value for a product species, or an estimated maximum 
for an intermediate species. 

Since 


k(T) 


AT 


e3q>(-T^/T) 


(6) 


the first bracketed term on the right side of equation (3) can be rewritten 


nco 

I ^(T) p(T) dT = 

f 

A(K,t+K) exp( J 

*'—00 ^ * 


L \ 1 2/J 


p(t) dt = k(t) 


(7) 
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where K- = and = T • • In similar fashion, the second bracketed 

I max min 2 min 

term on the right side of equation (3) can be written 


r r 





(8) 


where . = c . •• c . and 

j,max 3, min 


For computational purposes it 
the form 


^4,j 

is useful to rewrite the above relationships in 


w. 


= Zw. 


1,1 


1,1 




2 C- 

r A, A^ 


(9) 


where Z^, the temperature amplification ratio, is equation (7) divided by equa- 
Z^, the species amplification ratio, is equation (8) divided by the 

where c- is the mean value for species j; and 


tion (6); 
product 


^7V ^7V ' 

'"l ’^2 


z = Zj.Z^. Thus, 


represents the ratio of the turbulent mean value of the Arrhenius reaction-rate 
constant to its value calculated by inserting the mean temperature into equation (6), 
and Z^ is the ratio of the turbulent mean value of the product Ct, Cj, to the 


product 


^Ai"^A2 


""a ""a • 


Utilizing this formulation requires the selection of an appropriate pdf for the 
temperature (eq. (7)) and for the joint pdf of the species concentrations (eq. (8)). 
Ihe selection of the pdf is discussed in the following sections. 


Alternative Probability Density Functions for Modeling the Effect of 

the Temperature Fluctuations 

Ihree probability density functions for the temperature were discussed in refer- 
ence 17. These were the clipped Gaussian (ref. 12), the beta (ref. 18), and the ramp 
(ref. 19). The ramp probability density function has been observed in some regions 
of turbulent shear flows wherein temperature histories exhibit ramp-like behavior 
upon which high-frequency Gaussian fluctuations are superimposed (ref. 19). The 
principal conclusions reached in reference 17 concerning these three probability 
density functions are: 

- 2 

1. For all values of t and t* , the temperature amplification ratio Z^ is 
greater than unity; that is, the effect of the turbulent temperature fluctuations is 
to increase the effective Arrhenius reaction-rate constant. 

2. Values of Z^ obtained using the clipped Gaussian pdf are substantially the 
same as those obtained using the beta pdf. Differences between the values obtained 
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for develop principally at low and high values of the mean temperature t, which 

occur in regions where little reaction is occurring and where the reactions are near- 
ing completion, respectively. This generally close correspondence in the behavior of 
the temperature amplification ratio Z^ for the clipped Gaussian and beta pdf leads 
one to choose the beta pdf for practical calculations because of the relative ease 
with which the beta pdf can be employed in computer programs. 

3. Values of Z^ are substantially higher for the ramp pdf over the entire 
range of mean temperature t than they are for either the clipped Gaussian or beta 
pdf. 

An evaluation of the nonreactive flow-field data presented in references 20 
to 23, wherein ramp-like behavior of the tenperature fluctuations was observed in 
certain regions, leads to the following criterion for deciding whether to employ the 
beta pdf or the ramp pdf in equation (7) in a particular region of a turbulent flow 
field: 

- 2 

If 0.35 < t < 0.75 and t* > 0.08, then the temperature fluctua- 
tions are likely to exhibit ramp-like behavior. Elsewhere, the 
beta pdf should be used in equation (7). 

Another probability density function, not included in reference 17, is the 
"most-likely” pdf (ref. 5). Although it is computationally more difficult to utilize 
than the beta pdf, joint (multivariable) forms of the most-likely pdf can be written. 
As pointed out in reference 18, a physical basis exists for the use of the beta pdf 
to describe the turbulent fluctuations of many flow-field quantities. Hence, it was 
reasoned that if an essential equivalence could be demonstrated between the beta pdf 
and the most-likely pdf in its one-variable form, this would suggest the use of a 
bivariate form of the most-likely pdf to describe the species fluctuations. 

The most-likely pdf (ref. 5), in its one-variable form, is given by 

p(t) = CO) 


where the constant coefficients , and obtained from simultaneous 

solution of the following constraint equations: 

P p(t) dt = 1 

0 

P t p(t) dt = t 

0 

(t - t)^ p(t) dt = t'^ 

0 

Results obtained for using equation (10) in equation (7) are shown in 

figures 1 to 3. Parameter values for the cases cited in figures 1 to 3 are in 
table I. The most-likely pdf is compared with the beta pdf in the figures. Ihe two 
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pdf's yield quite similar results over the examined ranges of t and t* , and both 
exhibited similar sensitivities to the activation energy and temperature exponent. 
Again, this similarity favors the computationally simpler beta pdf in large-scale 
computer programs. However, as described above, the essential equivalence between 
the beta and most-likely pdf has important implications in selecting an appropriate 
joint pdf for the species concentrations in equation (8). This equivalence is exam- 
ined in the next section. 


Two-Variable Probability Density Function for Modeling the Effect 
of the Species-Concentration Fluctuations 

As noted in the preceding section, the similarity in the results obtained for 
the temperature amplification ratio using the experimentally verified beta pdf 

and the most-likely pdf suggests that a bivariate (or joint) form of the most-likely 
pdf serves to model the effects of the species-concentration fluctuations. To obtain 
values for the species amplification ratio in equation (9), the following bivar- 

iate most-likely pdf is used in equation (8): 


P 




exp(\o + 









( 11 ) 


where the value of q, called an a priori probability in reference 5, was, in accor- 
dance with reference 5, set equal to unity. The values of the constant coefficients 
^0' ^1' ^2' ^3 equation (11) are obtained from the simultaneous solution 

of the following constraint equations: 



The dimensionless mean species concentrations r and r are, as with 

— 2 . ^1 ^2 

t and t* , quantities which are calcu lated from transport equations describing the 

turbulent flow field. The covariance r* r* is determined from 

A, Aj 




r * r • = o * 


( 12 ) 


8 



where, as above, the dimensionless mean-square species-concentration fluctuations. 



and 



are calculated using appropriate transport equations 


(e. g. , ref. 


and where p* is the correlation coefficient, 
measure of the degree of linearity between r^ 


Ihe correlation coefficient 

and r^ • Values of p* 

^2 


P* 

near 


reflect a high degree of linearity, and values near zero indicate a lack of 
linearity. Positive values of p* indicate that as r^^ increases, r^^ 


increases; negative values of 
decreases. 


p* indicate that as r^^ increases. 


1 


2 ), 


is a 

±1 


From the data in reference 24, values between -0.9 and -1.0 were determined for 
the correlation coefficient p* between the tenperature and the hydrogen concentra- 
tion throughout the flow field. It was concluded that correlation coefficients of 
this magnitude appear to be realistic in at least some types of flows. Hence, for 
the present calculations, results were obtained for p* = ±0.9. 

“72 

Results for Z obtained using p* = -0.9 over ranges of r , r ' 

*1 ^2 

and r^ are shown in figure 4. In all cases < 1« 'Ihis implies that, in 

general, for reaction-rate terms involving species whose concentrations are changing 
in opposite directions (one increasing and one decreasing, which leads to p* < 0), 
the effect of the turbulent species fluctuations is to reduce the reaction rate. In 
the literature, this effect has often been referred to as unmixedness. Ihis unmixed- 
ness effect is seen in figure 4 to predominate (low values of the species amplifica- 
tion ratio Z ) when the mean concentrations of the species are low and turbulent 
fluctuations are large. 

On the other hand, figure 5 (where p* = 0.9) illustrates that for reaction-rate 
terms involving species whose concentrations are changing in the same direction 
(either both increasing or both decreasing, which leads to p* > 0), the effect of 
the turbulent species fluctuations is to increase the reaction rate through an 
enhanced mixing effect. The effect is again greatest at low values of the mean spe- 
cies concentrations and at high values of the turbulent species fluctuations. Ihese 
results suggest that these mixedness /unmixedness effects are likely to be greatest in 
the high -shear mixing regions of turbulent diffusion flames, where turbulent fluctua- 
tions are likely to be greatest. 


In writing equation (3), statistical independence was assumed between the tem- 
perature fluctuations and the species fluctuations. To examine this assumption, the 
results presented herein were conpared with the results obtained for an amplification 
ratio calculated assuming full statistical dependence between the temperature and 
species fluctuations. Ihe results for this latter amplification ratio Z^t,r^^,r^^^ 

were obtained by employing the following three -variable form of the most-likely pdf: 








3 A. 


X ±r 
4 A. 


X(-tr 
5 A^ 




(1 3) 
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where, as in equation (11), q, the a priori probability, is set equal to unity. 

values of \ , \ , ..., X are obtained from known values of t, r , r , 

^ T o A, A„ 


The 


1 


t'r' , t'r' , and r' r* by employing the following constraint equations: 
^1 *2 *1 ^2 


J’ /’ /’ ,r^ \ 

0 0 0 \ 1 2 / 


dt dr, dr, = 1 
"^1 ""2 


P f f t p/t,r^ , 

0 0 0 




dt dr, dr, = t 

Ai A^ 


f f f r P(t,r ,r \ 

0 0 0 1 V 1 2/ 


dt dr, dr, = r, 
A, Aj A, 


/' /' /' Pj P(‘-‘=A '"^A ) 

0 0 0 2 V 1 2j 


dt dr, dr, = r, 
""l ^2 ^2 


r f f (t - 1) 
0 0 0 


;A, - ^A,) P(-A -A,) 


dt dr^ dr = t'r* 
""l '^2 ''l 


f f (t -i)L \p 

0 0 0 \ 2 2 




dt dr_ dr_ = t*rl 
1 2 2 


i' r I' is - is - 


tS'S) 


dt dr dr = r' r! 

Ai A^ A^ A^ 


A typical result is presented in figure S, where ZCt^r, ) is the fully 

^2 

dependent amplification ratio calculated using equation (13), and Z = Z^Z is the 
result obtained assuming statistical independence between the temperature fluctua- 
tions and the species fluctuations. The results for this case are sufficiently 
similar in both magnitude and trend to support the assumption of statistical indepen- 
dence. In the next section, results are presented for a supersonic hydrogen-air 
diffusion flame for which experimental data were available. 


APPLICATION OF THE MODEL TO A SUPERSONIC, AXISYMMETRIC FLOW 


The chemical reaction model described in the previous sections was incorporated 
into the two-dimensional, parabolic computer program described in references 2 



and 25. This program solves parabolic partial differential equations for the trans- 
port of the following quantities: mass, momentum, energy, species concentrations, 

turbulent kinetic energy, dissipation rate of turbulent kinetic energy, species con- 
centration fluctuations, and temperature fluctuations. The program includes detailed 
chemical kinetics for describing chemical changes. The axisymmetric flow field exam- 
ined (shown in fig. 7) consists of a supersonic hydrogen jet injected into a super- 
sonic coaxial airstream. For these calculations, the eight-reaction H 2 -O 2 chemical 
kinetic mechanism is shown in table II. The values of the forward reaction-rate 
constants are also in table II. 

To utilize the portion of the model described herein dealing with the effect of 
the species -concentration fluctuations, a criterion is required for deciding whether 
to employ a positive or negative value for the correlation coefficient p*. That is, 
at a particular location in the flow field one must decide whether the concentrations 
of the species involved in the reaction are mutually increasing or decreasing 
(p* = 0.9) or whether one species is increasing in concentration while the other is 
decreasing (p* = -0.9). To ej^edite the calculations, the point in the flow field 
(at a given axial location) where the hydrogen and oxygen were present in stoichio- 
metric proportion was selected as the location for assigning particular trends to the 
species concentrations. For example, for radial values greater than that at the 
stoichiometric plane, water and the free radicals O, H, and OH were assumed to 
decrease in concentration as radius increased. Inside the stoichiometric location, 
these species were assumed to increase in concentration from the centerline to the 
stoichiometric point. On the other hand, H 2 was assumed to decrease monotonically 
and O 2 was assumed to increase monotonically from the centerline outward along the 
radius. As an example, for reaction 5 in table II, p* = 0.9 in both the forward 

and backward directions. For reaction 8, H 2 and O were assumed to move in opposite 

directions from the centerline to the stoichiometric location (p* = -0.9) and in the 
same direction (p* = 0.9) outside the stoichiometric location. However, for reac- 
tion 8 proceeding from right to left, the correlation coefficient is always positive. 

Results for the centerline hydrogen concentration calculated using instantaneous 
kinetics are compared with results obtained by incorporating the model proposed 
herein and with experimental values in figure 8. Slight improvement in the calcu- 
lated values using the turbulence reaction model may be seen, particularly in the far 
field. Radial pitot-pressure profiles at several axial locations are shown in fig- 
ure 9. Again, a slight improvement is brought about by the inclusion of the proposed 

^model in the far field in the region between the centerline and the stoichiometric 
location. Ihis improvement occurs as well for the radial species-concentration pro- 
files shown in figure 10. Ramp-like temperature fluctuations were not encountered in 

2 

any of the calculations described, because t* never exceeded 0.08. 

The results obtained from the flow-field calculations presented suggest a number 
of model refinements which are recommended for further investigation. Ihese refine- 
ments include the requirement for a more precise criterion to decide upon the sign of 
the correlation coefficient p*. For example, rather than utilize the stoichiometric 
H 2 -O 2 location as an arbitrary turnover point, the computer program could monitor the 
gradients of the species concentrations at previously calculated adjacent points to 
decide upon the sign of the correlation coefficient p*. Ultimately, the magnitudes 
of the concentration gradients could be used to assign values to p* between ±1, 
rather than using fixed values for p*. 



The criterion developed for observing ramp-like temperature fluctuations was 
based upon data obtained from nonreacting flows. Additional numerical experimenta- 
tion based upon alternative criteria for introducing temperature-ramp values of the 
temperature amplification ratio needs to be undertaken so that recommendations 

can be made for reacting flows. 

Additional calculations are required using chemical kinetic mechanisms other 
than the eight-reaction mechanism in table II. For example, reference 2 documents 
the improved correlation between calculated and experimental values obtained using a 
larger chemical mechanism for hydrogen-air combustion. Ihe effects of the model 
presented here on a larger-scale chemical mechanism must be assessed. 


CONCLUDING REMARKS 

It has been shown that a computationally viable, three-variable (temperature and 
species concentrations) model for reacting flows which accounts for the effects of 
turbulence on chemical reaction rates results in some improvement in the agreement 
between the theoretical prediction and the data from an axisymmetric, supersonic- 
combustion, hydrogen-fueled experiment. Ihe theoretical prediction of the flow 
parameters and species concentrations shows best agreement with the experiment in the 
field (greater than 20 fuel jet diameters (x/d.^^) downstream of the point of 
injection). Recommendations are made for what are"^considered necessary improvements 
to the model. These include improvements in methods for selecting a correlation 
coefficient, alternative criteria for introducing temperature -ramp values of an 
amplification factor, and improved correlations by additional comparisons with 
experimental data. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 28, 1 984 
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TABLE I.- TEST CASES 



TABLE II.- CHEMICAL KINETIC MECHANISM AND RATE CONSTANTS 


Reaction 

Forward rate constant^ 

A 

B 

E 

1. H2+M^H+H+M 

5.5 X lo’’® 

-1 .0 

51 087 

2. 02 + M;|:0 + 0 + M 

7.2 X 10 ''® 

-1.0 

59 340 

3. H 2 O + m;Joh + h + m 

5.2 X 10^^ 

-1.5 

59 386 

4. oh + m;|;o + h + m 

8.5 X 10 ^® 

-1.0 

50 830 

5. H 2 O + 0 t oh + OH 

5.8 X lo''® 

0 

9 059 

6. H 2 O + H ^ OH + H 2 

8.4 X 10 "'® 

0 

10 116 

7. O 2 + H t OH + 0 

2.2 X 10 '’^ 

0 

8 455 

8. H 2 + 0 t OH + H 

7.5 X 10 ''® 

0 

5 586 


a. B 

The form of the rate constant is k = AT exp(-E/RT) 

with k in cm /mole-s. 
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t = 0.5 


' — Beta and 
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t = 0.7 


(c) t = 0.5 and 0.7 


Figure 1.- Continued 













































(a) = 6.56. 

Figure 9.- Radial pitot-pressure profiles. 
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(b) = 1 3.80. 

Figure 9.- Continued. 
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(c) = 20.00. 

Figure 9.- Continued. 
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(d) x/dj^^ = 26.20. 
Figure 9.- Concluded. 


39 













1. Report No. 2. Government Accession No. 

NASA TM- 85746 

3. Recipient's Catalog No. 

4. Title and Subtitle 

A MODEL FOR REACTION RATES IN TURBULENT REACTING FLOWS 

5. Report Date 
May 1984 

6. Performing Organization Code 
505-43-83-03 

7. Author(s) 

Wallace Chinitz and John S, Evans 

8. Performing Organization Report No. 
L- 15679 

10. Work Unit No. 

9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665 

11. Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Memorandum 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 

14. Sponsoring Agency (^de 


15. Supplementary Notes 


Wallace Chinitz: The Cooper Union, New York, New York. 

John S. Evans: Langley Research Center, Hampton, Virginia. 


16. Abstract 


To account for the turbulent temperature and species-concentration fluctuations, a 
model is presented on the effects of chemical reaction rates in computer analyses of 
turbulent reacting flows. The model results in two parameters which multiply the 
terms in the reaction-rate equations. For these two parameters, graphs are presented 
as functions of the mean values and intensity of the turbulent fluctuations of the 
temperature and species concentrations. These graphs will facilitate incorporation 
of the model into existing computer programs which describe turbulent reacting 
flows. When the model was used in a two-dimensional parabolic-flow computer code to 
P^®*3ict the behavior of an experimental, supersonic hydrogen jet burning in air, some 
improvement in agreement with the experimental data was obtained in the far field in 
the region near the jet centerline. Recommendations are included for further 
improvement of the model and for additional comparisons with experimental data. 


page) 21. No. of Pages 

43 

For sale by the National Technical Information Service. Springfield. Virginia 22161 


22. Price 
AO 3 


18. Distribution Statement 

Unclassified - Unlimited 


Subject Category 07 


17. Key Words (Suggested by Author(s)) 

Parabolic computer program 
Supersonic combustion 
Turbulence modeling 
Turbulent reacting flows 

20. Security Ctassif. (of this 
Unclassified 


19. Security Qassif. (of this report) 
Unclassified 


NASA-Langley , 1984 























National Aeronautics and 
Space Administration 


THIRD-CLASS BULK RATE 


Washington, D.C. 
20546 

Official Business 

Penalty for Private Use, $300 


Postage and Fees Paid 
National Aeronautics and 
Space Administration 
NASA-451 





POSTMASTER: 


If Undeliverable (Section 1S8 
Postal Manual) Do Not Return 




